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The measurement of the Shapiro time delay in binary pulsar systems with highly inclined orbit 
can be affected both by the motion of the pulsar’s companion because of the finite time it takes a 
photon to cross the binary, and by the gravitational light bending if the orbit is sufficiently edge-on 
relative to the line of sight. Here we calculate the effect of retardation due to the companion’s 
motion on various time delays in pulsar binaries, including the Shaipro delay, the geometric lensing 
delay, and the lens-induced delays associated with the pulsar rotation. Our results can be applied 
to systems so highly inclined that near conjunction gravitational lensing of the pulsar radiation by 
the companion becomes important (the recently discovered double pulsar system J0737-3039 may 
exemplify such a system). To the leading order, the effect of retardation is to shift all the delay 
curves backward in time around the orbit conjunction, without affecting the shape and amplitude 
of the curves. The time shift is of order the photon orbit crossing time, and ranges from a second 
to a few minutes for the observed binary pulsar systems. In the double pulsar system J0737-3039, 
the motion of the companion may also affect the interpretation of the recent correlated interstellar 
scintillation measurements. Finally, we show that lensing sets an upper limit on the magnitude of 
the frame-dragging time delay caused by the companion’s spin, and makes this delay unobservable 
in stellar-mass binary pulsar systems. 

PACS numbers: 95.30.Sf - 95.85.Bh-97.60.Gb-97.60. Jd-97.60.Lf-97.80.-d 


I. INTRODUCTION 

Timing study of binary radio pulsars has provided 
some of the best tests of the general relativity (GR) to 
date. Among the effects that have been probed in differ¬ 
ent systems are gravitational redshift, precession of the 
periastron and orbital decay due to gravitational wave 
emission. Strong evidence for the geodetic precession - 
GR manifestation of the spin-orbital coupling of the bi¬ 
nary - has been found in several systems (e.g., 0 , 0 )- In 
binaries with highly inclined orbits, it has been possible 
to measure the Shapiro delay caused by the propagation 
of the pulsar s^nal in the companion’s gravitational field 
(e.g., d, 0, d). 

In systems that have nearly edge-on orbital orienta¬ 
tion, like the recently discovered double pulsar system 
J0737-3039 0 , one may also look for the effects related 
to the gravitational bending of the pulsar radio beam as 
it passes close to the companion around the moment of 
the pulsar’s superior conjunction ( 0 , 0 , 0 , 0 ). Light 
bending modifies the conventional Shapiro delay formula 
widely used in pulsar timing studies and introduces addi¬ 
tional geometric time delay ( 0 , 0 ). Bending also gives 
rise to a lens-rotational delay associated with the rotat¬ 
ing beam of the pulsar and induces distortion of the pulse 
profile shape near conjunction 0 . 
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Our previous studies of the lensing effects in binary 
pulsar systems ( 0 , 0 ) adopted a “static lens” approx¬ 
imation, in which the companion motion is neglected. 
Kopeikin & Schafer showed that the orbital motion 
of the pulsar companion gives rise to a Uc/c correction (vc 
is the companion velocity) to the conventional Shapiro 
delay formula and such correction might be observable 
under some circumstances. The work of [l?i| . however, 
neglected light bending and thus cannot be applied to 
the situation where the light ray passes the companion 
within a few Einstein radii. On the other hand, only 
for such highly inclined systems are the lensing effects 
potentially observable ( 0 , 0 ). 

In this paper we study the effect of companion mo¬ 
tion on the time delays associated with the gravitational 
light bending near the companion. These include both 
the lensing corrections to the delays independent of the 
pulsar spin - Shapiro and geometric delays (see gnu, and 
the lens-rotational time delays (gEl- Our formulae are 
general and valid in both the “strong lensing” and the 
“weak lensing” regimes. We also provide a simple and 
general interpretation of the influence of the lens motion 
on the behavior of various timing signals in binary pulsar 
systems (gvu which has not been previously discussed. 

Qualitatively, we expect that near binary conjunction 
the effect of moving lens amounts to a time shift in vari¬ 
ous time delay curves. It takes finite amount of time for 
the pulsar signal to travel across the binary system dur¬ 
ing which the position of its companion changes. Based 
on simple geometric considerations, this time shift is 
(Mp/M)a||/c, where Mp is the pulsar mass, M is the 
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total mass of the system, and a|| is a separation between 
binary components projected along the line of sight [see 
eq. 0 ], In the limit Mp 0 this retardation delay goes 
to zero since companion does not move. Our explict cal¬ 
culations and general expressions of various time delays 
confirm this simple expectation, see m 

For completeness, in m we discuss the effect of frame 
dragging due to the companion’s spin on the time delay. 
This t^e of delay was previously investigated by [llf . 
da, m, M, who disregarded light bending in their 
studies. We incorporate the lensing effect in our calcula¬ 
tion of the frame-dragging delay, and show that 

this delay is unobservable in the binary pulsar systems 
with stellar mass companions. 


II. LENSING BY MOVING PULSAR 
COMPANION 

Consider a binary system consisting of a pulsar with 
mass Mp and a companion with mass Me (the total mass 
is M = Mp + Me). Let Vp, be their position vectors 
with respect to the binary center-of-mass, and Rp, Rc be 
their projections on the plane of the sky. We also define 
r — Vp — Tc (the position vector of the pulsar relative to 
its companion), = —np-r (the projection of r along the 
line of sight; here np is the unit vector from the binary 
barycenter to the observer), and R = Rp —Rc = np x (r x 
np) (the projection of r in the sky plane). For a binary 
orbit characterized by the semimajor axis a, eccentricity 
e, longitude of periastron uj and inclination angle i, we 
have 


r = jrj = a(l — e^)/(l-I-ecosc^i), (1) 

r|| = rsinisin'i/), (2) 

i? = |R| = r(l — sin^ i sin^ '0)^^^, (3) 


terms of the binary orbital parameters we have 

&P = [rcosijj — d(sin'0 -I- esinw)]^ 

-I- [rsin^/^cosz -I- dcosi{cos'4’ + ecosw)]^ . (6) 


Here, to characterize the effect of lens motion, we have 
introduced the parameter d given by 


d = 


c-\/l — e 


Hbosinz 
'c(l -I- e sinw) ’ 


( 7 ) 


where Oc = a{Mp/M) is the companion’s semi-major axis 
and rib = {GM/a^y^^ is the binary angular frequency. 
Note that in eq. m and the remainder of the paper, 
all variables are evaluated at the time of pulse emission 
t = te- 

Because of light bending, the actual impact parameter 
of the ray (i.e. the minimum-approach distance between 
the radio beam and the companion), 6, differs from 6p. 
The bending angle, including the leadin g-or der correction 
due to the companion motion, is (e.g.. |15|. 


/ np-VcNdGMc 
“ “ c / c 26 ’ 


( 8 ) 


The Ve/c correction to the classical light-bending formula 
only slightly affects the value of the companion’s Einstein 
radius, and will be dropped in the remainder of the paper. 
We are only interested in lensing around the superior con¬ 
junction of the pulsar, during which 6p ~ i? ^ r|| ~ ay is 
well satisfied. Lensing gives rise to two images (positive 
and negative) of the source located at b± = 5±bp/6p in 
the plane of the sky with respect to the companion’s po¬ 
sition. From the lensing equation, b — bg = aa^\, we find 
the impact parameters of the two images: 


b± 


1 

2 




( 9 ) 


where tp is the true anomaly measured from the ascend¬ 
ing node of the pulsar and (p = tp — uj is the orbital true 
anomaly measured from the periastron. To describe grav¬ 
itational lensing of the pulsar signal by its companion, it 
is also useful to introduce 


0 || = a| sini|(l — e^)/(l-I-esinw), (4) 

which is the value of ry at the superior conjunction of the 
pulsar {tp = tt/2). 

A photon (radio pulse) emitted by the pulsar at time 
t = te from the position R(te) passes through the lens 
plane at time ti ~ tg -|- ry/c. In the absence of light 
bending, the impact parameter of the ray is 


bo 


Rp(te) - Rc{tl) = R{te) + Rc{te) “ Rc{tl) 


R{te) 


Ve{te) 

n—— 


no X 



X no 


(5) 


where Re is the Einstein radius given by 

Re = {2Rga\iy^y Rg = 2GMe/c^. (10) 

The image amplification A± = d{b'^)/d{bQ) is given by 
the standard expression 

+ 2 1 

^± = 7 ;— / 9 with u = bo/RE- (H) 

2uVu^ + 4 2 

Clearly, the motion of the companion affects the lens¬ 
ing images only through the “retardation” term in eq. ©, 
arising from the finite time ry /c that photons spend trav¬ 
eling between the pulsar and the lens plane of the com¬ 
panion. 


III. GEOMETRIC AND SHAPIRO DELAYS 


where Vc = dRc/dt = no x (vc x no) is the projection Lai & Rafikov 0 gave the expressions for the geometric 
of the companion velocity Vc in the plane of the sky. In delay and the lensing-corrected Shapiro delay for highly 
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inclined binary pulsar systems based on the static ap¬ 
proximation. Here we revise these expressions taking into 
account the motion of the companion. 

To order 0{vc/c), the expression for the geometric de¬ 
lay assumes the usual form (see Appendix A) 

(A<)geom = ^(^) , (12) 

where the image displacement is given by 


which agrees with the no lensing result obtained in M 
(see also [13). 

Not too far from the binary conjunction, when R ~ 
a|j[(AV')^ -I- (Az)2]i/2 < ay, eq. d reduces to 



+ ^ln[a(l-e2)] . (16) 


Ab± =b±-bo = ^ [±^bl + ml- 6o) . (13) 

The leading-order Vc/c correction to (At)geom enters 
through the expression for &o [see eq. ©]• Compared 
to the the static limit (3 the only difference is in using 
bo instead of R in eq. m- 

The easiest way to illustrate how the lens motion af¬ 
fects the Shapiro delay is to switch to a reference frame 
co-moving with the comDanion|39|. In this frame, the 
companion’s gravitational field is static and the expres¬ 
sion for the Shapiro delay derived in Q [see their eq. 
(5)] can be used. The expression involves ry and b. Be¬ 
cause of the aberration of light, the direction of pho¬ 
ton propagation in the comoving frame n' is given by 
n' ~ no — Vc/c. Thus, instead of ry one has to use 
rJi = —n'-r = ry -l-R- Vc/c. Also, the difference between 
b'j. (the impact parameter of the ray in the comoving 
frame) and b± is proportional to (uc/c)^. Thus, includ¬ 
ing the leading order Vc/c correction, thegravitational 
(Shapiro) delay is given by (cf. eq. (5) of [3)01 


(A<)„. = - (l - A 


X In 


Ra 


(ry + R • V,/c)2 + bl-^-R- V,/c 


-I —- In [a(l — 
c 


(14) 


Note that the prefactor (1 — ng • Vc/c), which was absent 
in the original formula of [3 , is not captured in our simple 
analysis given above. A formal derivation of eq. d is 
given in Appendix A. 

In the “weak lensing” limit, when bo ^ Re [or equiv¬ 
alently, [A'tpY + (Az)^ » (i?_E/oy)^ ~ (^/c)^, where 
Aij) = — and Az = z —7r/2], we have for the positive 

imagejiil, - Re/^o = 2ay('u/c)(i?B/&o) < a||(z;/c), 
where v ~ (GM/ay)^/^. Thus, b+ ~ R— ryVc/c. Keep¬ 
ing only terms up to order 0{vc/c), the numerator of 
the expression inside the logarithm in eq. d reduces to 
r — ry — R • Vc/c = r -|- r • [no — ng x (vc x no)/c]. As a 
result, in the absence of lensing, eq. d becomes 


(At)grav = - (l - 


Rn 


X In 


Ra 


r -I- no • r — 


c / c 

no X Vc 


(no X r) 


-I —- In [a(l — e^)l . 
c 


( 15 ) 


where / ~ 1 — R • Vc/(cry) ~ 1. Note that the regions 
of validity of the limiting cases described by eqs. itTBl) 
and d have a significant interval of overlap [(w/c)^ <C 
(AV>)2 + (Ai)2<l]. 

IV. LENS-ROTATIONAL DELAYS 

Because of light bending of the pulsar signal by the 
companion, the direction of emission n at the pulsar po¬ 
sition differs from ng, with the deflection vector (An )e = 
n — no given by 

Ab± /A6±\ no X [(r-ay Vc/c) X no] 

(An)L =-= ^-, 

ay \ bo J ay 

(17) 

where Ab± = b± — bo = Ab±{ho/bo) [see eqs. @ and 
dj- Again, we are interested in the lensing effect only 
around the superior conjunction of the pulsar, so that 
ry ~ oy. Since pulsar signals are due to the beamed 
emission of the neutron star (as opposed to radial pulsa¬ 
tion), such a change in the emission direction gives rise 
to time delays and pulse profile distortion which depend 
on the pulsar spin period ([3)E3' 

In the following we will use the description of the pul¬ 
sar spin geometry from |18| (see their Fig. 1) in which 
the pulsar spin axis Sp (the unit vector along the pulsar 
spin direction) is specified by two angles: C.p is the angle 
between Sp and the light-of-sight vector ng, and rip is the 
angle between the ascending node of the orbit and the 
projection of Sp on the sky plane. 


A. Longitudinal Time Delay 

Variation of the emission direction (Anj^ causes a shift 
in the equatorial longitude $ of the emission direction in 
the corotating frame of pulsar (counted in the direction 
of rotation): 


^ (An)L ■ (sp X ng) 
|Sp X ng |2 


(18) 


This longitudinal shift corresponds to a change of the 
emission phase, leading to a longitudinal time delay 
(At)i = A$/r2p (Up is the angular frequency of the pul¬ 
sar; the time delay is positive for signal arriving later). 
This delay is independent of the specific emission pattern 
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(as longas it is rigidly rotating around Sp) and is given 
by (cf. 13) 


(At)L 


/ (r-a||Vc/c) ■ (sp x rip) 

\ bo J f^poylsp X nop 


(19) 


J 


The longitudinal time delay shifts the pulse profile uni¬ 
formly in time without affecting its shape. Evaluating 
eq. (PI explicitly using the pulsar spin geometry de¬ 
scribed above, we find 


(At)L = 




sin rjp cos ip — cos i cos r]p sin ip 
rip sin ^p 



sin? 7 p(sin '0 + esinw) + cos i cos rjp{cos ip + e cos uj) 
rip sin (pp 


( 20 ) 


where bo is given by eq. ©■ Setting d = 0 recovers 
the expression for the longitudinal lens-rotational delay 
obtained in [3 under static approximation. 


B. Pulse Profile Variation and Latitudinal Time 
Delay 


Variation of the pulse emission direction (An)^ also 
results in the change of colatitude (p of the emission vector: 


(AC)l = - 



Sp ■ (An)^ 

|sp X no| 

no X (r - ayVc/c) 


(sp X "o) 

ISpXnoP ^ ^ 


ing vector model (RVM, see ^3) commonly adopted for 
inferring the pulsar spin-magnetic orientation from the 
pulse profile and polarization data. In the course of pul¬ 
sar rotation around Sp, a given emission cone is traversed 
by the line of sight to the observer no only if the cone’s 
half-opening angle p satisfies p > |Cp — a|. Crossing the 
edges of the cone by no leads to the leading and trailing 
emission episodes observed at Earth. 

In the RVM, the latitudinal shift (AOl [due to no ^ 
no -t- (An)i] results in a variation of the full width 2$o 
of the pulse profile (corresponding to a given emission 
cone), with |3 


A$o = - 


(AQl 

sin (pp tan xo 


( 22 ) 


In general, any shift of the (co)latitude of the emission 
vector leads to a pulse profile variation because the pulsar 
radio emission is thought to be tied to a specific geometric 
pattern on a celestial sphere rigidly rotating around Sp. 
Since the efficient particle acceleration leading to radio 
emission is thought to take place at the pulsar magnetic 
polar caps, the radio emission pattern is expected to be 
connected with the position of the pulsar magnetic axis 
m. 

Here, for illustrative purposes, we assume that the ra¬ 
dio emission pattern forms a set of coaxial cones with cir¬ 
cular cross sections and symmetry axes coincident with 
the magnetic axis m, which is inclined at an angle a 
with respect to the pulsar spin axis Sp (see Eig. 3 of |3 
for details of the adopted pulsar spin-magnetic geome¬ 
try). This emission pattern forms the basis of the rotat- 


Here xo is the angle (on the celestial sphere) between 
the arc connecting np and Sp and the arc connectingno 
and m at the edges of the pulse (see Eigs. 3 & 4 of j3). 
In the framework of the RVM, xo is directly related to 
the position angle of the linear polarization of the pulsar 
emission, and is given by [23 


tan Xo 


sin a sin $o 

cos a sin Cp — cos $o sin a cos (pp ’ 


(23) 


The pulse contraction/expansion A$o directly trans¬ 
lates to a latitudinal delay of the arrival time of various 
pulse components. This delay is (At)^'^*^ = A$o/Hp for 
the leading component of the pulse profile fi^j- Combin¬ 
ing eqs. m and 122) we find 


J 


X 



cos r]p cos Ip + cos i sin pp sin ip 
flp sin (pp tan xo 



(A^)(iat) ^ (AQl ^ / Ab± \ 

^ HpSinCptanxo \ bo ) 

cos i sin rjp (cos ip + e cos uj) — cos rjp (sin ip + e sin uj) 
ftp sin ^p tan xo 


(24) 


Setting c? = 0 in eq. (| 22 | reproduces equation (23) of [ 3 , 


which was derived in the static approximation. 
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Since depends through xo on phases $o of dif¬ 

ferent pulse components, the latitudinal lensing time de¬ 
lay leads to inhomogeneous dilation/contraction of the 
pulse profile around the superior conjunction of the 
pulsar 12 . This is in contrast to the longitudinal shift 
of the emission vector, which only causes homogeneous 
displacement of the whole pulse profile in time without 
introducing any distortions of its shape. This difference 
of behaviors led to propose timing of the individual 
features of the pulse profile as a means of detecting the 
latitudinal shift. Such a procedure yields a larger infor¬ 
mation content than the standard timing of the whole 
pulse profile (assumed to have fixed shape) which can 
only reveal the longitudinal time delay. 


V. FRAME-DRAGGING TIME DELAY 

Frame-dragging delay (At)FD owes its existence to the 
gravitomagnetic effect of the companion spin angular mo¬ 
mentum Sc. In order of magnitude, the metric perturba¬ 
tion due to Sc is ~ Sc/r'^ (in units where G = c = 1), the 
additional time delay for a ray with impact parameter b 
is then (At)FD Sc/b. Given the small magnitude of 
this delay, in its subsequent consideration we will neglect 
the effect of the companion motion and will use static ap¬ 
proximation accounting only for lensing effects. We have 
(see e.g., 0) 


(At)FD 



dVph * 


Sc X {rph - rc) 
\TCph - TcP 


(25) 


where Vph is the position vector for the photon trajectory 
(see also Appendix El). Assuming that the observer is 
at infinite distance from the system and evaluating the 
integral in ll^ we find 


(At)FD 


2G no • (Sc X R) 
b±R 



.(26) 


As (At)FD may only be noticeable in nearly edge-on sys¬ 
tems and at the conjunction, one can use b± rn. We 
specify the orientation of the companion’s spin vector Sc 
by the angles C,c and rjc (analogous to angles (jp and rjp 
used to describe the orientation of Sp) and its magnitude 
by the dimensionless spin parameter d < 1 according to 



(27) 


where we have introduced a fiducial unit of time 


Rl 

(At)““ = 

cKe 

« 1.44 X 10“^^s 



(29) 


Note that (At)FD never diverges, not even for Az = 0. If 
one were to neglect the light deflection by the companion 
(^has been done in all previous studies, see p^. 0, 
[T^ . Q), (At)FD would still be given by equation ll^ 
but with R instead of b±. This would lead to a divergent 
(At)FD as i —> 7r/2 and ip ^/2. In the limit R ^ 
Re, eq. (1^ reduces to the corresponding result of 0 
obtained neglecting lensing. 

From eq. (1^ . it is clear that (At)p“ is the maximum 
possible value of the frame-dragging time delay since 
6+ > Re- This fundamental limit may be approached 
if the companion is a maximally spinning black hole with 
a certain spin orientation and i = 7r/2, ip = 7r/2. Note 
that for the negative image, b- < Re and (At)FD exceeds 
(At)““ roughly by RE/b-, (neglecting various geomet¬ 
ric factors and assuming d = 1), but its magnification 
A_ = (6_/i?F;)^/[l — {b_/RE)'^] rapidly goes to zero as 
the time delay increases. Thus, observations of the neg¬ 
ative image are not feasible in practice. 


VI. APPLICATION TO HIGHLY-INCLINED 
SYSTEMS 

The effects of gravitational lensing and companion mo¬ 
tion discussed in previous sections are of interest only for 
nearly edge-on binaries, those with |Az| = |z — 7r/2| 1. 

In particular, for the “strong-lensing” systems (those 
with |Az| ^ Re/qw), the image shift at the orbital con¬ 
junction {ip = 7r/2), Ab ~ Re ~ a,\\v/c, is of the same 
order as the apparent change of the photon impact pa¬ 
rameter due to the companion motion, |bo —R| = a||14/c 
[see eq. ©]■ Therefore it is important to account for both 
lensing and companion motion properly for such systems. 


A. Effect of Companion Motion 

For a highly inclined binary system (|Az <C 1|) around 
the orbital conjunction {\Aip = \ip — 7r/2| <C 1), the ef¬ 
fect of lens motion on the various time delays studied in 
previous sections has a simple interpretation, as we now 
show. Expanding eq. (0 for |Az|, \Aip\ <^1, we find 


We then find 


fRE\ 


U± j 

-l-(Ai)^ [l -1- (d/a||)ecosa;]^| 


sin r]c cos ip — cos rjc sin ip cos i 
(1 — sin^ Ip sin^ z)i/^ 


1/2 


(28) 


I [Alp (d/a||)(l -I- esinw)]^ -I- (Az)^| (30) 
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Comparing with the static limit, = i? ~ 

a|| [(AV')^ + we see that 6olv> - &o‘'"*l</'+Aj/'ret) 

where we have introduced a phase shift 


A^/iret = —(1 + esinw) 

Oil 


flbac (1 + e sinw) 
c Vl - 


(31) 


Clearly, in the presence of the lens motion, the photon 
impact parameter is minimal not at Aip = 0 (exact con¬ 
junction) but at —A'0ret) which occurs slightly before the 
conjunction. This is a direct manifestation of the retar¬ 
dation effect due to the finite time it takes for the photons 
emitted by the pulsar to propagate through the binary, 
during which the lensing companion’s position changes 
(see 0, P). 

As mentioned in m near the conjunction the Shapiro 
delay (A<)grav is described by eq. itT^ . It is then obvious 
from 03 and that lens motion affects (Atjgrav and 
(At)geom only through its effect on bo which enters b± in 
these formulae. Thus, 

(At)grav |-!/> — (At)gjav li/'+Ai/iret) 

(Atjgeomli/) — (At)geomlV’+Ai/’reti 

for |Ai|, \Aif\ < 1, where and (At)ggo*ja are 

the Shapiro and geometric delays calculated in the static 
limit i. 

For the spin-dependent time delays near conjunction, 
one finds using eqs. (HOI) and 


(At)L - 


m 


cosr]pAi — s'mrip [Atf + ((i/a||)(l -I- esinw)] 


flp sin 
Ab± 


(Ai)dat) _ _ 


sin? 7 pAi -|- cos rjp [Aip + ((i/a||)(l -I- esinw)] 


rip sin (p tan xo 


,(33) 


.(34) 


It is again clear that 


(Atjil^ ~ ( At)|)’"% + AV.ret > 

(At)^'"')|^ ~ (35) 


where (At)|)®'* and are the longitudinal and 

latitudinal light-bending delays calculated in the static 
approximation (i). 

Equations and show that all delays discussed 
here are related to their static analogs via {At){'ip) ~ 
(Ai)stat(^ -H AV'ret) for |Ai|, |AV>| < 1. Thus, the ef¬ 
fect of the companion motion is to simply shift the delay 
curves (At)®*”^‘(■)/)) homogeneously in phase by —Aif^et 
without affecting the amplitude or shape of these curves 
(to v/c accuracy). This phase shift is equivalent to mov¬ 
ing the time delay curves calculated in the static approx¬ 
imation earlier in time by 


Ay^ret (1 - 
rib (1-I-esina;)^ 


Mp a 1 — e^ 
M cl + e sin u; 


(36) 


One can see that as advertised in this time shift is 
indeed equal to (Mp/Mjay/c, i.e., the light crossing time 
of the binary at conjunction weighted by the mass ratio 
Mp/M. 

To illustrate our results, consider the double pulsar 
J0737-3039 system, consisting of the millisecond pulsar A 
with spin period Pa = 22.7 ms and mass Mp = I.337M0, 
and the normal pulsar B with mass = 1.25Mq play¬ 
ing the role of the lensing companion 0). This 

nearly edge-on system has Pb = 0.1023 d, e = 0.0878, 
and a = 8.784 x 10^° cm. The phase shift is Aip-cet « 
1.2 X 10“^, and time shift Atj-et ~ 1-4 s. In Fig. 1 we dis¬ 
play the timing delays discussed in this paper for J0737- 
3039 - the combined geometric and Shapiro delays given 
by eqs. ca and m, and the rotational lensing delays 
given by eqs. and lEU ~ together with the image 
magnification factor given by eq. lllll , which depends on 
bo and is thus also affected by the lens motion. We choose 
i = 90.29° ([i^) for this figure, which allows lensing ef¬ 
fects to manifest themselves (see a, a). To evaluate 
the spin-dependent delays (At)i and (At)^^*^*^ we adopt 
as an example C, = 50°, rj = 45° consistent with the po¬ 
larization measurement of |24| . Compared to the static 
results, all these curves are shifted back in time from the 
moment of conjunction. 

Although the example in Fig. 1 refers to the “strong 
lensing” case (|Ai| ^ v/c), we emphasize that for any 
highly-inclined system (|Ai| <C 1), the effect of the com¬ 
panion motion can be captured adequately by the “re¬ 
tardation shift” as described above. In Table I we list 
Ay^i-et and Af^et for binary pulsar systems with measured 
Shapiro delays. One can see that the value of Af^et ranges 
from several seconds in compact systems to several min¬ 
utes in wide binaries. Detection of the retardation ef¬ 
fect seems easier for the latter but would ultimately de¬ 
pend on the accuracy of timing observations of a given 
svstem|4^. 


B. Effect of Frame-Dragging 

Previous conclusions on the observability of the frame¬ 
dragging time delay have been over-optimistic as they 
have always been obtained neglecting the light bending 
effect. For example, for a binary pulsar system with 
Me = 10 Mq, 0=1, and a|| « 9.4 Rq (correspond¬ 
ing to binary period of 1 d), Wex & Kopeikin m found 
that (AfjpD reaches 1 fis for \i — 7r/2| « 0.05°, while 
our eq. (1^ gives (At)p“ « 0.15 fis. This difference 
is easily understood by noticing that for this set of pa¬ 
rameters, Re « 2 X 10^ km is considerably larger than 
the assumed unperturbed impact parameter of the radio 
beam, a|||i — 7r/2| « 5.7 x 10^ km. Thus, light bend¬ 
ing cannot be neglected in the calculation of (Atjpo for 
this system. For binary pulsars with neutron star com¬ 
panions, (At)pD is completely negligible (contrary to the 
claims in Refs. la and M)- 

Another problem with stellar-mass binaries is that 
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FIG. 1: The amplification (top panel), the combined geometric and Shapiro delay (the second panel), the rotational-lensing 
delay (the third panel), and the latitudinal lensing delay (the bottom panel) of the dominant (“plus”) image of the pulsar A 
signal as a function of the orbital phase in the double pulsar PSR J0737-3039 system. The longitude is measured from the 
superior conjunction of pulsar A (when A is exactly behind pulsar B). The inclination angle is chosen to be 90.29°. For the 
third and bottom panels, the angles ^ = 50°, rj = 45° and tanxo = 0.08 are used. In each panel, the solid curve corresponds 
to the static limit (see Fig. 1 of Q), while the dashed curve includes the effect of the companion motion. 



















TABLE I: Binary pulsars with measured Shapiro delays. 


PSR 

ms 

Pi’’, d 

a/c, s 

Me, Mq 

sin 2 

A'0ret 

Atret 7 S 

Ref. 

B1534-H2° 

37.9 

0.42 

7.6 

1.345 

0.975 

9.2 X lO"'* 

8.87 

[2^ 

40737-3039° 

22.7 

0.1 

2.93 

1.25 

0.999987° 

11.8 X lO"'* 

1.39 

[3j 

B1855-H09'^ 

5.4 

12.3 

60 

0.23 

0.9995 

3 X lO”'* 

50.8 

m 

J0437-4715'' 

5.8 

5.7 

37.6 

0.24 

0.679 

4.2 X lO"'* 

32.6 

[27] 

J1713-H0747‘* 

4.57 

67.8 

192 

0.28 

0.9 

1.7 X lO"'* 

158 


J1640-h2224‘* 

3.16 

175.5 

369 

0.25 

0.99 

1.3 X lO""* 

313 

m 

41909-3744'* 

2.95 

1.53 

15.3 

0.2 

0.998 

6.4 X lO"'* 

13.4 

[^ 

40751-H807'* 

3.4 

0.26 

5.8 

0.19 

0.83 

14.9 X lO-* 

5.32 

[29] 


“Pulsar spin period, Pp = 27rjQp. 

^Orbital period of the binary, = 27r/n^. 

"^Pulsar companion is a neutron star. 

^Pulsar companion is a white dwarf. 

^Value of sini from scintill atio n measurements of |^. Shapiro 
delay suggests sini = 0.9997 13(1 . 

(At)FD is degenerate with the longitudinal lens- 
rotational delay (At)/,. This was first pointed out in 
M under the weak light bending approximation^^, but 
it is easy to see that this degeneracy holds even in the 
case of strong lensing; Ab± satisfies the lens equation 
Ab± — R%/b±, and the ^dependent factors in eqs. ll^ 
and are analogous]^ (the only difference is in rjp 
used in the former and rjc in the latter). As a result, these 
two delays are covariant and in practice it is not possible 
to separate their contributions unless there exists some 
additional information on their relative magnitudes such 
as the knowledge of the pulsar spin-magnetic orientation 
from polarization measurements. 

These examples lead us to conclude that measuring 
frame-dragging effect in the stellar mass pulsar binaries 
is impossible if the timing precision is to remain at the 
level of 1 /iS. In principle, shrinking binary (reducing 
a||) boosts up (At)p“, but this comes at the cost of 
reducing the binary lifetime due to gravitational wave 
emission El. If a pulsar is found around a high-mass 
black hole (such as the one in the Galactic Center, see 
|3l|l in a highly inclined orbit, the chance of measuring 
frame-dragging delay may be better: in this case (At)p“ 
not only is large but also strongly exceeds (At)^. 


VII. DISCUSSION 

In this paper we have studied the combined effects of 
gravitational lensing and companion motion on various 
time delays in binary pulsar systems. Our study im¬ 
proves upon previous work on the lensingeffect by go¬ 
ing beyond the static approximation (@,0)- We show 
that for highly inclined systems, the companion motion 
affects the Shapiro delay and other lensing-related de¬ 
lays by shifting the delay curves backward in time by the 
amount ~ Oc/c (where Oc is the semimajor axis of the 
companion). 


To detect this effect on time delays, one should look 
for a systematic displacement of the peak of the Shapiro 
delay curve from the exact moment of the superior con¬ 
junction of the pulsar. To this end, one must be able 
to pinpoint the exact moment of the conjunction and 
to nail down the peak of the Shapiro delay curve. This 
may not be so easy since in practice one bins the tim¬ 
ing data within rather large time intervals to get a good 
measurement of the Shapiro delay. Thus, one would need 
repeated measurements of many binary pulsar orbits to 
be able to reach the accuracy needed for measuring the 
retardation effect. Retardation effect on lensing delays 
may be detectable in the double pulsar system J0737- 
3039, but the feasibility of this measurement depends 
on the currently uncertain inclination of the system and 
is compromised by the eclipse of the millisecond pulsar 
near its superior conjunction by the magnetosphere of its 
companion (i, M)- 

In the case of the double pulsar J0737-3039 system, 
the retardation effect due to the companion motion not 
only affects timing of the millisecond pulsar but may also 
influence the interpretation of the scintillation measure¬ 
ments of this system near conjunction |23|. Indeed, the 
correlated measurements of the interstellar scintillations 
of the two pulsars provide information about the sys¬ 
tem orientation and properties of the interstellar turbu¬ 
lence only when the positions of both pulsar projected 
on the sky plane are accurately known at each moment 
of time. In the static approximation the projected sepa¬ 
ration of pulsars is R, while it becomes bp [see eq. 0] 
when the retardation effect is taken into account. For 
the currently inferred inclination angle from the scintil¬ 
lation measurement itself (|i —90°| = 0.29°±0.14°, [^j. 
both the retardation effect and the light bending distor¬ 
tion of the apparent path of the lensed millisecond pulsar 
@ have similar orders of magnitudes and contribute at 
the level of ~ Re to the projected pulsar separation at 
conjunction. As for this system Re = 2550 km is not 
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very different from the inferred minimum projected ap¬ 
proach distance of the two pulsars, « 4000 ± 2000 km, 
one expects that both effects may at some level affect 
the interpretation of the correlated scintillation measure¬ 
ments of 1^. Whether this can explain the discrepancy 
between the values of the system’s inclination obtained 
using the Shapiro delay {i = 88.7° ± 0.9°; 0, [s^) and 
scintillations is not clear at present. 
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A Lorentz transformation shows that the correction to 
eq. (IA8II is of order and will be neglected. The 

photon trajectory rph(t) satisfies = 0, which gives 

dt = [l- 2$(1 - 2vc • k)] dl, (A9) 

where drph = d/ k (with k the unit vector of photon 
propagation) and dl = (drph • drph)^/^. This corresponds 
to an effective index of refraction 1 -|- (2Mc /1 x' |) (1—2vc • k) 

(e.g., HI)- 

Consider a photon emitted from the source (the pulsar) 
at time tf, (near the binary conjunction), and arriving to 
the observer at to- The impact parameter of the photon 
isb = bo + Ab [see eq. 03] ■ Integration of eq. 1A9II along 
the ray path yields 

to-te = |rph(to)-rph(4)| + (At) geom -b(At)grav, (AlO) 


APPENDIX A: FORMAL DERIVATION OF THE 
SHAPIRO DELAY FORMULA FOR MOVING 
LENS 


where the 1st term gives the usual Roemer delay, the 2nd 
term is the geometric delay given by equation m, and 
the 3rd term is the Shapiro delay represented by 


The metric produced by the companion (mass Me) in 
its rest frame (“primed” frame) assumes the standard 
form (we set G = c = 1 in this appendix) 

ds^ = -(1 + 24>)(df')^ + (1 + 24>)-idx' • dx', (Al) 

where <I> = — Mc/|x'|, and x' is the 3-vector measured 
from Me- In the “lab” frame (“unprimed”) comoving 
with the barycenter of the binary, where the lens moves 
with velocity Vc, the metric can be obtained by a Lorentz 
boost. To order 0{vc), it is given by 

ds^ = —(l-|-2<I>)dt^-|-8$Vc-dxdt-|-(I —2$)dx-dx. (A2) 

The above result can also be obtained using the linearized 
theory of gravity (e.g., [ISl): The metric coefficients pro¬ 
duced by a source with density p{t, x) and velocity v(t, x) 
are 


(At)grav = ydZ ^^(I - 2vc • k), (All) 

where the integration is along the ray path and x' = 
rph(^') - J'c(i') - rph(f) - Veit). 

Without loss of generality, let us assume that the pho¬ 
ton passes through the lens plane at t = ti = 0. Then 
the photon position vector is rph(0 = rph(O) -I- ctk, 
with k = no for t > 0 and k = n = ng -I- An, 
where ng is the unit vector pointing toward the ob¬ 
server and An = Ab/a||. The lens position vector is 
Te = rc(0)-|-Vct = rc(0)-|-Vct-l-U||tno, where uy = ng-v^, 
and Vc is the projection of Vc on the sky plane. Thus, 
along the ray x' = (b — Vet -I- crcAnt) + {c — U||)tno, 
with (7 = 1 for t < 0 and cr = 0 for t > 0, and 
b = b± = rph(O) — rc(0). Taking the integral in (lAllll 
we have to order v'^ 


where e 


goo — + 2C/ -I- 0{e^), 

(A3) 

go, = -4U,+G(e5/2), 

(A4) 

Qij = (1 -l- 2U) + 0{e^), 

(A5) 

0{vl, Mc/r), and 


t;= /■'’(' 1 ,/"■='') dV, 

J |x-x'| 

(A6) 

./ X-x' 

(A7) 


Taylor expansions show that the retardation effect enters 
the potentials o nly in the order 0{Ve), and eq. is 

recovered (e.g., [l7|l. 

Let the position vector of the companion (lens) be 
rc(t), and that of the ray be rph(t). Along the ray, the 
potential $ = — Mc/|x'| should be evaluated at 


x' = rph(i') - i’c(i') ^ rph(t) - Veit). (A8) 


Afi = 2Mc(l -I- U||) In -I- \/s'^ + 6^^ 
+2iWc(l -l- U||) 

X In -I- \/s'^ + 52 ^ 


i=So(l-ii|| )-b.Vc 
i=-b.Ve 


—b-(Vc—cAn) 


s=Se(l—'i^ll)—b-(Vc —cAn) 


(A12) 


where Sg = etc, So = cto- Since Sq —> 00 , the upper limit 
only amounts to a constant which can be dropped. Thus 
to leading order 


Ati — —2il(fc(l Uj|) 

X In fs -I- y/s^ -I- b'^'] (A13) 

V / s=Se(l— f II )-b.(Ve—cAn) 

Note that Sg = ctg is determined by ng • 
[rp(fg) -h (0 - 4)cn] = ng • rc(0) = ng • [rc(tg) - Vgtg], 
which gives Sg = ctg = —(l-t-uy) ry, where ry = —ng-r(tg) 
and r = Tp — rg. Thus the Shapiro delay (Af)grav = 
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(1 — 2v||)Aii is given by 


(Ai)grav = -2Mc(l - no • Vc) 

X In ^(ry + b • Vc — cb • An)^ + 

-r|| - b • Vc + b • An] , (A14) 


where all quantities are evaluated aX t = te- 

There are some differences between eq. (IXni) and 
eq. (I14II (apart from the trivial additive constant term 


in the latter). These differences appear in the expression 
inside the logarithm, and are of order (uc/c)^. This is 
expected since our calculation includes only the leading- 
order Vc/c correction to the “static” Shapiro delay for¬ 
mula, while higher-order terms are neglected. Indeed, it 
is easy to show that eq. imi reduces to eqs. m or 
(fTHIl in the appropriate limits. We conclude that both 
eqs. m and iTmil give the same Vc/c corrections to the 
Shapiro delay, and they should in practice give the same 
quantitative result for (At)gi.av 
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